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Abstract. We study the kinetics of nonlinear irreversible fragmentation. Here 
fragmentation is induced by interactions/collisions between pairs of particles, and 
modelled by general classes of interaction kernels, and for several types of breakage 
models. We construct initial value and scaling solutions of the fragmentation 
equations, and apply the " non- vanishing mass flux "criterion for the occurrence of 
shattering transitions. These properties enable us to determine the phase diagram 
for the occurrence of shattering states and of scaling states in the phase space of 
model parameters. 
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Fragmentation is a phenomenon of breakup of particles into a range of smaller 
size particles. It is naturally found in a wide variety of physical systems, ranging 
from comminution, breakup of grains, bubbles, droplets, polymer degradation, 
disintegration of atomic nuclei, etc. Fragmentation may occur through external forces, 
spontaneously, or through interactions/collisions between particles. The subject has 
been widely studied [T]-[l0]. 

We are mainly interested in collision-induced nonlinear fragmentation as caused 
by binary interactions. Such systems can be described by the time evolution of c(x, t), 
which is the number of particles of mass or size a; at a given time t, or alternatively 
by its moments M n (t) = J °° dxx n c(x, t). Quantities with similar properties appear in 
coagulation processes. In either case the total mass is conserved, M[t) = M\{t) = 1, 
while the total number of particles, N(t) — Mo(t) — J dxc(x,t) is not. In irreversible 
coagulation, the mean particle mass, s(t) = M/N(t), increases monotonically, and 
may lead to a finite time singularity at t = t c , the gelation transition, characterized 
by the appearance of an infinite cluster containing a finite fraction of the total mass, 
A(t) = 1 — M(t) (order parameter), where A(t) ^ for t > t c . Alternatively, the 
gelation transition is characterized by a non-vanishing mass flux A(t) = —M(t) from 
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finite size particles (sol) to the infinite cluster (gel) [H]-[I3], i.e. a violation of mass 
conservation. 

In irreversible fragmentation the reversed scenario occurs. Here s(t) is 
monotonically decreasing, while the overall mass is conserved. In these systems a 
finite time singularity may occur at t c , the shattering transition. It is characterized by 
a non- vanishing mass flux, A(t), i.e. the rate at which massive particles are converted 
into mass-less infinitesimals or fractal dust [U EJ [7] . If A(t c ) is finite this transition 
has the character of a continuous phase transition, described by the order parameter 
A(i) = 1 — M(t), as in gelation [TJ]- In case A(i c ) = oo admass instantly 'evaporates' 
from the system; the transition is called explosive and also referred to as a first order 
transition [8]. 

Smoluchowski's coagulation - fragmentation equation [llj gives the basic mean 
field description for reversible and irreversible coagulation (2j E[ [13] and fragmentation 
processes [T]- [TU] in terms of the time evolution of c(x.t) in spatially uniform 
(well stirred) systems. In irreversible fragmentation or coagulation the system is 
described by a nonlinear coagulation rate, in combination with a spontaneous linear 
fragmentation rate and/or a collision- or reaction-induced nonlinear fragmentation 
rate. The system does not reach a steady state, but at asymptotically large times 
the distribution function c(x 1 t) approaches under rather general conditions to the 
standard scaling form, which describes the typical £— dependence around the mean 
particle size s(t), which is steadily decreasing. 

The occurrence of shattering has been addressed only partially in the case of 
collision-induced nonlinear fragmentation. It shows a behavior, qualitatively different 
from spontaneous (linear) fragmentation. Furthermore, the special cases analyzed so 
far are not necessarily generic, but appear to be borderline cases. In this letter we 
study the occurrence of shattering for general classes of fragmentation models within 
the framework of the nonlinear fragmentation equation and we analyze its peculiarities 
and point out the parallels with gelation. 

Collision-induced irreversible fragmentation can be described at the mean field 
level by the nonlinear fragmentation equation with a collision term / composed of a 
loss and a gain term g], 

dc(x,t)/dt = I(x\c) = -c{x)f™dyK(x,y)c(y) 

+ I™ d y Jo° dzb ( x \y) K (y> z ) c (y) c ( z )- W 

Here b(x\y) is a conditional probability, describing the distribution of outgoing 
fragments of mass x, given that a particle of mass y breaks [TJ [2j O 0]. One 
distinguishes: (i) deterministic or splitting models [HIE], where a particle breaks into 
two equal fragments, hence b(x\y) = 25(x — y/2), and (ii) stochastic models, where 
a fragment of random mass x breaks off from a particle of mass y. As mass is 
conserved in a single breakup event, the outgoing fragment distribution has to obey 
the homogeneity requirement, b(x\y) = y~ 1 b(jj). For simplicity, we take the standard 
form b(s) = ((3 + 2)s< 3 [M], obeying, 

J^dxxb(x\y)=y , N = J* dxb(x\y) = |±f. (2) 

For physical reasons the mean number of outgoing fragments satisfies N > 2 which 
implies — 1 < (3 < 0. Binary breakup corresponds to {3 = 0. In Eq.JT]) we consider 
binary interactions, where the kernel K(x, y) describes the interaction rate of pairs of 
particles (x,y). It may further contain a factor p(x,y), which gives the probability 
that breakage indeed occurs, and may depend on the masses {x,y). If p is constant, 
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it can be absorbed in the time scale. Of importance in our analysis is also the rate 
equation for the cumulative mass, M(x, t) = J dyyc(y,t), as can be derived from 
Eq(l). It reads, 

M{x, t) = -/~dy/ °° dz yB(Z)K(y, z)c{y)c{z) (3) 

with B(s) = Jq duub(u) — s@ +2 . In applications of the nonlinear Smoluchowsky 
equation a variety of collision kernels K has been proposed for processes induced 
by interacting particles [TTJ [T3]. Because of mathematical simplicity, kernels of 
the sum-product form: K(x,y) — l,x p + y p , (xy) p , x p y q + x q y p , etc have been 
extensively studied in coagulation processes. Physically motivated kernels are e.g. 
K = (x~ a + y~ a ")(x" + y 13 ) with a = j3 = 1/d for interaction rates among diffusing 
particles, or K ~ (R x + R y )"'~ x for ballistic collision rates in a <i-dimensional system. 
Here R x ~ x x l d is the radius of a particle of mass x, and K a geometrical cross-section. 
In most cases of physical interest, especially at limiting particle masses (x <C s(t) or 
x 3> s[t)), the kernels are continuous and homogeneous, i.e. K[ax, ay) = a x K(x, y) = 
a x K(x,y), and b[ax\ay) — a x b(x\y) with A' = — 1. Kernels for coagulation can be 
classified by two exponents [13] . i.e. K(x,y) ~ x p y q if x <C y, where p < q and 
X = p + q, with p > (class I), p — (class II) and p < (class III), with the physical 
restrictions A < 2, q < 1 |13) . This (p, q)- classification appears to be relevant for 
nonlinear fragmentation as well, as we will show. 

The breakage probability, p{x,y), defines three different types of models 
depending on whether particle x or y breaks: (i) symmetric breakage, where a randomly 
chosen particle of the interacting pair (x, y) breaks [U[S], an d where p(x, y) — 1; (ii)i- 
breakage, where the larger particle breaks, hence p{x,y) = 9{x — y); (iii) S-breakage, 
where the smaller particle breaks and p{x,y) = 9(y — x); 9{x) stands for the unit 
step function. The corresponding nonlinear fragmentation equation for L-breakage 
is obtained from Eq.JT]) by replacing I(x\c) with Il(x\c), with K(x,y) replaced by 
Kl(x, y) — K(x, y)9(x~y), and similarly for S-breakage. Subsequently, we will discuss 
the nonlinear stochastic fragmentation equation for kernels of class I, II and III for 
symmetric, L- and S- breakage models. Regarding exact solutions of the nonlinear 
fragmentation equation very little is known, and mostly restricted to monodisperse 
initial conditions. The essential references are [U [8], where the former contains a 
representative list of the older literature. Ref.[4] analyzes the deterministic L- and S- 
breakage model, KL(x,y) — x p 9{x — y) and KL{x,y) — x p 9{y — x), and Ref.[Sj does so 
for both the deterministic and the stochastic breakage models with Kl(x, y) = 9{x — y) 
and K L (x,y) = 9(y - x). 

Regarding the structure of the nonlinear integral-differential equations (1) and (3) 
for cases where K{x,y) = a{x)a{y) is a general product kernel, it has been observed 
[1], that the nonlinear fragmentation equation can be transformed into a linear one 
with a new time variable r(t) that is related to the physical time t in a nonlinear 
manner. The functional form of r(t) determines whether a shattering transition is 
present or absent. So, to explain this dependence it is paramount to discuss how the 
initial solutions c(x, t) of the nonlinear fragmentation equation for a given c{x, 0) can 
be constructed from the initial solutions c(x, r) of the linear fragmentation equation. 
To this end, we analyze linear fragmentation, 

dc(x, t)/dt = -a(x)c(x) + f™dyb{x\y)a(y)c(y), (4) 

where a(x)c(x) represents the spontaneous or externally-induced linear breakup rate. 
Exact solutions c(x,t) are known for algebraic fragmentation rates, a(x) = x a for 
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all real a, and mono-disperse initial conditions, c(x,0) = 8(x — xq) [TJ [51 [31 [5]. 
These c(x,t)'s are the causal Green functions of Eq.Q with a monomer source 
6(t)S(x — x ) [5 \. So, c(x,t) = for all x > x at t > 0. In the sequel we set 
Xo = 1. Spontaneously fragmenting systems [15] with a > are non-shattering, i.e. 
the total mass M(t) = 1 at all times, and the total number of particles, Mo(i) < oo 
for all t < oo. Moreover, moments with n + [3 + 1 > exist, and evolve for large t as 
M n {t) ~ ft 1 "™)/" while those with n + (3 + 1 < are divergent. On the other hand, 
spontaneously fragmenting systems with a < 0, are shattering, and mass loss starts 
at the initial time. So, shattering occurs at t = t c = 0, where A(t c ) is finite; hence 
the transition is continuous. The possibility of an explosive shattering transition with 
A(t c ) — oo is never realized. All initial solutions, which have by definition t > 0, 
are non-mass conserving post-shattering solutions with t > t c = 0. They behave 
for small x as c(x,t) ~ A(t)x~ e with 6 — a + 2. Consequently, moments M n (t) 
with n < 1 + a are divergent for all times, and those with n > 1 + a decay for 
long times as M n {t) oc Ait) cx e^ t t^ +2 ^ a , including M(t). We also point out that 
for a > the exact solutions converge asymptotically to a standard scaling form, 
c(x,t) = (1/ s 2 (t))(p(x/ s(i))pQ-[9], where the scaling limit is formally defined as the 
coupled limit, t — > oo and x — * with x/s(t) kept constant. In this limit, where 
s ~ 1/Mo ~ t -1 /" — > 0,(i.e. the total number of particles Mq(£) diverges) the exact 
solution becomes, s 2 c(su, t) ~ (^(u) ~ w' 3 exp[— u"]. Those with a < do not approach 
a scaling form. Inspection of the exact solution as a — > at fixed (x, t) shows that 
c(x, t) = for all ir > x = 1, and reads for x < 1 (see Ref. [2]), 



where s(t) = e - *, and is the modified Bessel function of integer order n = 1 . 

This expression shows that the borderline case, a = 0, is exceptional, i.e. non-scaling 
and non-shattering. 

Let us now consider the nonlinear fragmentation equation for symmetric breakage 
with product kernel K(x, y) = a(x)a(y) — (xy) p and < A = 2p < 2. In this 
case Eq.([T]) is a quasi-linear equation for which exact initial value solutions can be 
obtained. It can be mapped onto the linear fragmentation equation with a — > p and 
t — > r, defined through dr = M p (t)dt. Consequently the mass distribution, c(x,t) 

and moments M n (r) = dxx n c(x,r) for mono-disperse initial values are known 
explicitly, and only r(t) needs to be determined in order to have the complete solution 
as a function of t. For p = a > 0, where K is a class I kernel, total mass is conserved for 
all r, and the moments for n ^= 1 behave at large r as M„(r) ~ T^ 1 ~ n ^ p . Furthermore, 
to have r as a function of t we need to invert the relation t = J Q ds/M p (s). If 
A = 2f> > 1, then t(r) ~ T 2_1 / p is monotonically increasing, the relation is invertible, 
and t — ► co as t — > oo. Hence, Mi(t) = 1 for all t, and there is no shattering and 
no divergence of Mo(t) at any finite time. However, if A = 2p < 1, then as r — > oo, 
t — > t c = J dr/Mpij) < oo. Consequently there exists a finite time singularity, 
t - (t c - t)-p/(i-2p) as t _> tc an d 

mass remains conserved only for t < t c and 
vanishes instantaneously at t c , where A(t c ) = — oo. At the same time all moments, 
behaving as M n (t) ~ (t c — t)( 1_ ™)/( 1_2 P) j either diverge or vanish. These are the 
hallmarks of an explosive shattering transition at t c where all massive particles are 
converted instantaneously into fractal dust. For A = 2p < (class III) the kernel 
K = (xy) p can be mapped on the linear equation through dr — M_\ p \(t)dt. Its 
moments M„(r) with n < 1 + a — 1 — \p\ do not exist. Consequently r(t) is not 



c(x,t) ^e- t {2t/\n(l/x)) 1 / 2 h 2(2t\n(l/x)) 1/2 



(5) 
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defined, and c{x, t) does not exist for mono-disperse initial conditions with class-Ill 
kernels. The same applies to scaling solutions. The corresponding class II kernel with 
p = a = {) 01 K= \ represents an exceptional point, as discussed below Eq.(5). The 
solution at p = (K = 1, class II) exists, is shattering, and c(x, r) is identical to the 
non-generic, non-scaling solution of the linear Eq.Q at a = 0. 

To determine possible scaling solutions we substitute the scaling ansatz c(x, t) — 
(l/s 2 )<p(u = x/s) in ((3J, and take the derivative, yielding the scaling equation for 
symmetric, L- and S-breakage (A = 0, L, S), 

{cp(u)/u?) ' = -(cp(u)/^ +1 )J^dvK A (u, v)<p(v), (6) 

where 7 is an arbitrary positive separation constant, ip(u) has to satisfy the boundary 
condition, u 2 ip(u) — as u — > 00, and Kq = K. For all three types of breakage models 
the evolution equation for the mean particle size is the same, ss~ x = —7. Its solution 
is, 

C (i +i) 1/(1 ~ A) forA > l,t -> 00 
s(t) ~ < exp[-7*] forA = l,t — y 00 (7) 

[ {t c -t) 1/{1 - X) forA < l,t <t c 

where to,t c — y 00 as A — ► 0. The appearance of the finite time singularity at t c 
indicates that shattering only occurs for A < 1. Systems with A > 1 are non-shattering 
[4]. Note that the scaling limit in the pre-shattering critical region is defined as the 
coupled limit: t ] t c and x — > with x/s(t) — constant . 

For sum-product kernels in symmetric breakage the rhs in Eq.® reduces to sums 
of powers u s , multiplied by coefficients m n — J Q duu n ip(u), which can be determined 
self-consistently. Specifically, for K(x,y) = (xy) p (p > 0, class I), one obtains from 
the rhs of (6) ip'( u ) = il){u)m p u p ^ 1 /j, where ij)(u) = ^p{u)/u^. This can be readily 
integrated to yield, 

(p(u) = Cvr exp[~u p m p /jp], (8) 

where C and 7 are determined by imposing normalization (mo = 1) and mass 
conservation (mi = 1). In order to get simpler analytic expressions, we use the 
invariance property that the scaled distribution, <p(u), obtained under the similarity 
transformation <p(u) — SQ 2 ip(u/ so) for an arbitrary constant so also satisfies Eq.© 
with fh\ = 1. This property allows us to fix 7 by setting m p /^p = 1, which is a 
self-consistency equation. With this choice the moments read m n — (C/p)T(b n ) = 
r(6„)/r(6i) where T(x) is the Gamma function and b n = (1 + (3 + n)/p. Hence, 
the scaling distribution function can be expressed as <p(u) = pu^e~ u /T((/3 + 2)/p). 
Solutions with a different normalization, e.g. mo = 1, are easily derived using the 
invariance property under similarity transformations. 

Similarly one derives that the size distribution for sum-kernels, K — x p + y p (class 
II) with p > 0, has the form (p(u) = Cu^er uP = pu^-^^h^" /T((f3 + 3)/2p) with 
[3 — (3 — m p /j. Here we impose the self consistent equation m /pj = 1 to determine 
the separation constant 7. Imposing mass conservation leads to the second equality 
in the previous equation. The moments of the distribution can then be computed; in 
particular m p — bojp where now b n = (1 + j3 + n)/p which implies (3 = ^((3—1). 
The exact scaling solutions for the symmetric breakage kernels (p,p) and (0,p) above 
approach different limiting forms as p — > 0. So, the analysis starting below Eq. §5§ 
shows that the if —kernel with (p, q) = (0, 0) is quite singular. In a similar manner 
the scaling solutions <p{u) for the geometric collision cross-section, K ~ (x 1 / 3 + y 1 ^ 3 ) 2 , 
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and closely related kernels can also be found, as well as the asymptotics of <p(u) for 
general class I and II kernels in all breakage models of type A = (0, L, S). In L- and 
S-breakage models for generic K no exact initial value or scaling solutions are known, 
except for the non-generic borderline case K = 1 in Ref.[5], which lacks standard 
scaling in the variable u = x/s(t) in all three breakage models. 

To analyze from a broader perspective the occurrence of shattering, we will 
focus on the behavior of the cumulative mass flux for vanishingly small masses. If 
lim^^o M(x, t) = M(t) is vanishing at t c , the system is non-shattering; otherwise 
there is shattering. If — oo < M{t c ) < 0, the phase transition is continuous, and 
c(x,t) exists for t > t c . If M(t c ) = — oo, the phase transition is explosive (first order), 
and c(x,t) does not exist for t > t c . 

For general K(x,y) in class I, II and III the mass flux M(x,t) can only have 
a non-vanishing limit for x —>■ if c(x, t) is of power law type, because the rhs 
contains the factor x^ +2 with — 1 < (3 < 0. So we propose the post-shattering ansatz 
c(x,t) ~ A(t)x- (x -> 0), and determine 6 such that M(t) ^ (see Refs.QU [13]). 
Moreover 9 < 2 because the total mass should remain finite. 

The evolution equation for fragmentation with S-breakage is described by Eq. ([3J 
with K replaced by K$- Inserting the ansatz above yields for small x, 

M(x, t) ~ -A 2 k{9)(x 3+x - 29 )/(29 + /3 — 1 — X) 

k{9) = dsK(l,s)s- e {l+p<9<2) (9) 

where K(l,s) ~ s p for s 3> 1. In case 6 — |(3 + A) the above small- a; limit 
yields a finite result for the mass flux M(t) = — A 2 (t)k(9) / (2 + (3), i.e. it allows 
the existence of a continuous shattering transition with a post-shattering solution of 
algebraic form for t > t c ; 9 < 2 implies A < 1. At the (unknown) shattering time 
t c mass conservation breaks down, and for t > t c there exists a non-vanishing order 
parameter A(t) = 1 - M{t) > with A(t) — A 2 (t). Eq. © includes also the special 
result, obtained in [8] for the S-breakage model with K = 1 and j3 = 0. 

It is remarkable that the post-gelation distribution, c(x,t) ~ A(t)x~ s , occurring 
in Smoluchowski's coagulation equation for A > 1, has the same exponent 9 = 
(3 + A)/2 P31 [12] as in the fragmentation process above; hence the close analogy 
between gelation and continuous shattering. Note that the value of exponent [3 has 
no influence on the existence of shattering. 

A similar analysis can be performed for symmetric and L-breakage models. In 
doing so we introduce a lower cut off ey on the ^-integral in Eq. (|3|), and take lim e ^o 
at the end of the calculations. Due to the physical restrictions on the allowed values 
for p and q, shattering is always explosive rather than continuous. 

From the properties discussed in this letter we can construct the phase diagram 
for symmetric breakage in the (p, q)— plane. It is restricted to the triangular region, 
spanned by (0, 0), (0, 1), (1, 0) and includes the boundaries. The region with < A < 1 
represent shattering systems, and the region with 1 < A < 2 represents non-shattering 
ones. The whole triangular region shows standard scaling in the variable u — x/s(t), 
except in the singular corner (0, 0). Regarding the phase diagram for L- or S-breakage 
the location of the left boundary (separatrix between "non-existence" and "existence 
of scaling solutions" ), including the singular point (0, 0), is unknown, and the behavior 
on it may be different from its right and left limit. From Ref.[8] it is known that a new 
type of scaling in the variable x/m*(t) appears at (0,0). Here m* is a characteristic 
mass, that cannot be defined a priori, but follows from a clever mapping of the 
fragmentation equation on the nonlinear equation for travelling fronts. Contrary to 
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models with symmetric breakage, which are quasi-linear, the scaling equations for 
L- and S-breakage are genuinely nonlinear, i.e. f"{u) = F(u,ip',(p), and the only 
solutions known are those for the singular point (0, 0). 

We have discussed the generic behavior of collision-induced irreversible 
fragmenting systems at the mean-field level. We have shown that the scenarios 
for nonlinear fragmentation are qualitatively different from those of spontaneous 
linear fragmentation. The behavior of the shattering transition depends both on the 
kind of fragmentation kernel and on the type of breakage. For symmetric and L- 
breakage, where the kernel K has a degree of homogeneity A < 1, shattering is always 
explosive, while S-breakage models show a continuous shattering transition, analogous 
to gelation. The existence of a transition does not depend on the details of the fragment 
distribution, b(x\y), i.e. on [3. Shattering in collision-induced fragmentation always 
takes place at a finite time t c ^ 0, as opposed to linear fragmentation where shattering 
occurs at t c — for a < 0. Contrary to gelation [T3], in class 111 kernels with symmetric 
breakage neither initial, nor scaling solutions exist. The solutions for fragmentation 
models with a fragment distribution, b(s) = for s < sq, b(s) ^ for so < s < 1, have 
scaling solutions if{u) exhibiting a log- normal distributions at small u [4]. 
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